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ABSTRACT 

In the present contribution, I summarize a systematic study of ISO and Spitzer 
mid-IR spectra of Galactic regions and star forming galaxies. This study quanti- 
fies the relative variations of the main aromatic features inside spatially resolved 
objects as well as among the integrated spectra of 50 objects. Our analysis im- 
plies that the properties of the PAHs are remarkably universal throughout our 
sample and at different spatial scales. In addition, the relative variations of the 
band ratios, as large as one order of magnitude, are mainly controled by the 
fraction of ionized PAHs. In particular, I show that we can rule out both the 
modification of the PAH size distribution and the mid-IR extinction, as an ex- 
planation of these variations. High values of the I6.2/I11.3 ratio are found to be 
associated with the far-UV illuminated surface of PDRs, at the scale of an inter- 
stellar cloud, and associated with star formation activity, at the scale of a galaxy. 
Using a few well-studied Galactic regions, we provide an empirical relation be- 
tween the I6.2/I11.3 ratio and the ionization/recombination ratio Go/ue -s/T^ss- 
Finally, I show that these trends are consistent with the detailed modeling of the 
PAH emission within photodissociation regions, taking into account the radiative 
transfer, the stochastic heating and the charge exchange between gas and dust. 

Subject headings: astrochemistry — (ISM:) dust — ISM: lines and bands — ISM: 
molecules — infrared: ISM, galaxies 



1. INTRODUCTION 

Polycyclic Aromatic Hydrocarbons (PAHs) are large molecules containing ~10 to 1000 
carbon atoms. They are commonly believed to be the carriers of the ubi quitous broad mid 



infrared features, centered around 3.3, 6.2, 7.7, 8.6 and 11.3 /im (e.g. lAllamandola et al. 
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19991 ). Due to their small size, these molecules are predomantly excited by single UV-photon 
events. At solar metallicity, ^ 15% of the tota l infrared (IR) luminosity is radiated through 



the aromatic features (e.g. IZubko et al.l 2004 ). On th e contrary, PAHs are underabundant 



in low-met alicity environments (iGalliano et al.ll2008al . for a complete study) 



It has long been surmised that the mid-infrared emission features provide a clear sig- 
nature of the interaction of far-UV photons with cloud surfaces and hence a probe of th e 
importance of massive star formation in a region (IGenzel et al.l Il998l : iPeeters et al.ll2004l ). 
The present work extends this by studying methodically the mid-infrared spectra of a wide 
variety of sources at different spatial scales, and investigating the physical origin of the PAH 
band ratio variations. 
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Fig. 1. — Qualitative demonstration of the aromatic feature variations in galaxies IjGalliano et al\\2008H) . 
Both panels show mid-infrared spectra normalised by the monochromatic Rux in the 11.3 fim feature. Left 
panel: integrated spectra of three galaxies; a blue compact dwarf, UM 448; a spiral M 51; a ULIRG 
NGC 4945. Right panel: spatially resolved spectra in three different regions in the starburst galaxy M 82; 
(A) is in the halo; (B) is at an end of the disc; (C) is in the starburst region. 
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SYSTEMATIC ANALYSIS OF THE OBSERVED AROMATIC 

FEATURES 

2.1. Sample and Mid-Infrared Spectral Decomposition 



We have constructed a large sample of mid-infrared spectra of Galactic regions (Hll 
regions, photodissociation regions, planetary nebulae), Magellanic Hll regions, and nearby 
dwarf, spiral and starburst galaxies. These sources were observed with one of the three 
following instruments: ISO/CAM, ISO/SWS, Spitzer/IRS. Several of these objects have 
been spectrally mapped (IC 342, M 17, M 51, M 82, M 83, 3 Doradus, and the O rion 
bar). Our sou r ces are the merging of the samples presented by iMadden et al.l (120061 ) and 
Galliano et al.l (|2008aJ). They present a wide range of properties (metallicity, star formation 
rates), as well as of spatial scales (resolution of ~ 0.1 pc in M 17 and the Orion bar, and of 
~ 0.1 kpc in external galaxies). Fig. [1] illustrates the diversity of these properties. 
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Fig. 2. — Demonstration of the two methods used to decompose all our mid-IR spectra, on the spectrum 



of M 82 iGalhano et al\\2008a) . Left paneh the continuum (yeUow Riled) is the sum of two modihed black- 
bodies; the aromatic features (blue shaded) are Lorentz profiles; the ionic lines (pink shaded) are gaussian 
functions. Right panel: the continuum (yellow Riled) is a spline function; the plateaus under the bands 
(empty areas) are also spline functions; the aromatic bands and ionic lines (blue shaded) are the residuals. 
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Quantifying the spectral variations demonstrated in Fig. [T] is not straightforward. In- 
deed, a significant fraction of the power radiated through the aromatic features comes from 
their wings. Moreover the accurate p r ofile of the bands is unce rtain, and has been m odeled 
differently by iBoulanger et al.l (119981 ) , iLi fc Draind (120011 ) and IVermeij et al. In or- 

der to test the robustness of our results, we have systematically analysed ou r spectra using 
two dif ferent methods (Fig. [2]). These methods are described in detail by iGalliano et al. 
( l2008bl ). It is important to note that the Lorentzian method (left panel of Fig. [2]) fits Lorentz 
profiles to the aromatic features, taking into account the flux in their wings. This method 
is feature-biased, since part of the small grain continuum can artificially be accounted by 
the wings of the bands. On the contrary, the Spline method (right panel of Fig. [2]), which 
integrates only the tip of the aromatic bands, is continuum-biased. 

Although the numerical values of the intensity of a given feature, measured with the two 
methods are different, our study shows that the order between several intensities is identical 
with the two methods. It indicates that the respective biases of the methods do not affect 
the general trends between band intensities. Therefore, in the following of thi s review, I will 
present only the results obtained with the first method and refer the reader to iGalliano et al. 
(l2008bl ) for the complete results. 



2.2. Correlations Between Aromatic Feature Intensity Ratios 

Fig. [3] shows the correlations between several aromatic feature intensity ratios derived 
from integrated spectra of galaxies and Galactic regions (left panels), as well as from spatially 
resolved spectra within M 82 (right panels). These figures clearly show that all types of 
sources and all types of apertures follow the same universal trends. The 6.2, 7.7 and 8.6 fim 
features appear to be tied together, while the ratios between these features and the 11.3 fim 
band vary by one order of magnitude. The relations involving the 8.6 /xm feature exhibit 
a larger dispersion, as this band is the weakest of the four considered here, its intensity is 
significantly affected by the silicate extinction feature, and it is blended w ith the powerfu l 



7.7 fim feature. A preliminary version of these correlations was presented by lGalliand (120041 ). 



In general, high values of I7.7/I11.3 are associated with intense star forming regions. Fig.|4] 
shows images of several spectral components of the edge-on Galactic photodissociation region 
(PDR) M 17. On each image, the Hll region is located toward the top-left corner, and the 
molecular cloud toward the bottom-right corner. The edge-on view of this region allows us to 
see the different strata of species, as demonstrated by the top left panel of Fig. |H This panel 
shows that the Ne ill ions are distributed on a ridge located closer to the H ll region than the 
Ne II ions. The top right panel shows that the Ne ill and the very small grain intensity have 
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Fig. 3. — Correlations between select aromatic feature ratios llGalliano et ali \2008ni) . l\ refers to the 
intensity of the band centered around A /xm. Left panels: the symbols correspond to the band ratios derived 
from spectra integrated over an entire source; the faint symbols are fits that we consider uncertain, due either 
to a low signal-to-noise ratio or to a weak feature-to-continuum ratio. Right panels: the symbols correspond 
to the band ratios derived from spectra of individual pixels within M 82; the three regions (A, B, C) are the 
ones shown in Fig. H In all panels, the grey stripe is the linear correlation to the observations ±1 x cr. 
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Fig. 4. — Spatial distribution of the various mid-IR components within the Galactic PDR M 17 
Cesarskv et al. IQOd: Galliano et al. 20081). Each image compares two spectral components; the inten- 



sity increases from dark red to light yellow for the color images; the contour levels are 5, 10, 20, 30, 50, 
70 and 90% of the peak intensity. The displayed components are the intensities of the [Neii]i2.8ipm and 
[Neiii]i5.56pm lines, of the continuum (Icont, integrated between 10 and 16 fim), of all the PAH bands IpAHj 
and the I7.7/I11.3 ratio. The shaded circle in the bottom right corner of each hgure is the beam size. 
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similar spatial distributions, while the bottom left panel shows that the PAHs are found 
more deeply inside the cloud, since they undergo photolysis/thermolysis close to the Hll 
region. Finally, the bottom right panel shows that the I7.7/I11.3 ratio peaks in the region 
where the [Ne Ill]i5.56;im and the continuum intensity Icont are maximum. Consequently, at 
the scale of an interstellar cloud, the I7.7/I11.3 ratio is maximum at the surface of the far-UV 
illuminated surface. 
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Fig. 5. — Spatial distribution of the aromatic features in M 82 (left panel) and M 51 (right panel), from 
Galliano et al. I \2008]\) . The color images show the intensity of all PAH bands (increasing from dark red to 
light yellow). The contours show the values of the I6.2/I11.3 ratio. The shaded circle is the beam size. 



Fig. [5] compares the spatial distributions of the intensity radiated by the aromatic bands 
to the I6.2/I11.3 ratio (following the same trends than I7.7/I11.3), within two external galaxies, 
M 82 and M 51. Contrary to Fig. HI the clouds are not resolved in this case. These images 
indicate that the I6.2/I11.3 ratio is enhanced where the PAH intensity is high, along the spiral 
arms and in the starburst regions. Therefore, at the scale of a galaxy, the I6.2/I11.3 ratio is 
maximum in massive star forming regions. 



3. ORIGIN OF THE AROMATIC FEATURE VARIATIONS IN GALAXIES 

3.1. Inventory of the Physical Processes Potentially Affecting the Band Ratios 

PAHs are stochastically heated, hence their thermal emission spectrum is determined by 
their temperature distribution. Fig. [6] demonstrates the effects of changing this temperature 
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Fig. 6. — Demonstration of the effects on the aromatic bands of changing the PAH temperature distribution. 
Each panel shows the theoretical spectrum of a collection of PAHs of various sizes N(j (number of carbon 
atoms; different colors). The total spectrum integrated over the size distribution is the black thick line. 
Top-left paneh we adopt the Zubko et al. 1 200-^) size distribution, and the Galactic diffuse interstellar 



radiation held (ISRF: \Mathis et al\\l98dl) : this case serves as reference for the other panels. Top-right panel: 
the size distribution is truncated; we exclude grains smaller than Nq ~ 1000. Bottom-left panel: the ISRF is 
harder (instantaneous burst of Myr). Bottom-right panel: the ISRF intensity is XISRF = 10^ times higher. 
The absorption efficiencies are a mixture of 50% neutral and 50% ionized PAHs, from \Draine & 
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distribution by varying different parameters. Fig. [6^ corresponds to the Galactic diffuse ISM. 
Fig. [Gt shows that the bands around 6.2, 7.7 and 8.6 fim decrease relative to the 11.3 fim 
band, if the small PAHs (smaller than Nq = 1000 carbon atoms) disappear. Indeed, small 
PAHs ffuctuate up to higher temperatures than larger ones. Hence, their short wavelength 
features emit more intensely. This case corresponds to environments where the small PAHs 
would photosublimate preferentially. Fig. [6t shows that by increasing the hardness of the 
interstellar radiation field (ISRF) or equivalently by increasing the average stellar photon 
energy, PAHs fluctuate up to higher temperatures, therefore boosting the short wavelength 
features. Finally, Fig. [6li shows that increasing the radiation density decreases the time that 
a given PAH spend on average at cool temperatures. However, this effect plays a role only 
for ISRFs which are Xisrf ^ 10^ times more intense than the Galactic one. 
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Fig. 7. — Demonstration of the effects on the aromatic bands of the extinction and of the charge of 
the PAHs. Both panels show the theoretical spectra with the \Zubko et al. honi) size distribution, the 
Draine & L\ 12001 ) absorption efficiencies, irradiated by the Mathis et al. \l 98 j) ISRF. Left panel: the 
reference spectrum (blue) is compared to the same spectrum cxtincted by a screen with an Ay ~ 10 (red). 
Right panel: the red and blue spectra correspond to neutral and ionized PAHs. 

Fig. [7] demonstrates additional effects controling the aromatic feature ratios. First, mid- 
IR extinction is dominated by a broad silicate feature centered around 9.7 yum. This feature 
absorbs significantly the 8.6 /im and 11.3 fim bands, but has little effect on the 6.2 and 
7.7 fim bands, as shown on the left panel of Fig. [3 Second, the 6 to 9 fim features are 
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Fig. 8. — Comparison between observed and theoretical trends of aromatic feature intensity ratios. Left 
panels: observed correlations within M 82 (right panels of Fig. the grey stripes are line hts ±1 x a to the 
observations. Right panels: presentation of the results of the modeling shown in Figs. [60 In each panel, 
the blue empty polygon with the circled corners shows the effect of varying both the minimum cut-off PAH 
size, -/V™™, and the ISRF intensity, XISRF, with the Galactic ISRF (labeled Galaxy). The values between 
parenthesis, close to each corners are the values of these two parameters at the corner. Second, the red empty 
polygon is the analog of the previous one, except that the ISRF is now a harder one (labeled Burst). Finally, 
the two filled polygons are the analogs of the two previous ones, except that we applied a screen extinction 
with Ay ~ 10. The various arrows show the sense of variation of the band ratios with each parameter. 
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enhanced for charged PAHs (right panel of Fig. [7]), since they are attributed to C-C mo des, 
while the 11.3 fim is attributed to peripheral C-H modes (e.g. lAUamandola et al.lll989l ). 



3.2. The Dominant Role of the Fraction of Ionized PAHs 

Fig. [H] compares the trend between band ratios observed within M 82 (similar trends 
are followed by the other sources; Fig. [3]) to theoretical trends. The left panels of Fig. [8] 
are another way to look at Fig. [3l Statistically, the observed I7.7/I6.2 and I8.6/I6.2 do not 
show significant variations and are uncorrelated with the I7.7/I11.3. The right panels of Fig. [8] 
represents the variations of the aromatic features from the models presented in Figs. 0- 
[71 What this figure shows is that the modification of the PAH temperature distribution 
(Sect. 13.11) is not consistent with our observations, since it would anticorrelate I7.7/I6.2 and 
I8.6/I6.2 with I7.7/I11.3. On the other hand, the extinction produces the right trend between 
I7.7/I6.2 and 17,7/111.3, but it anticorrelates I8.6/I6.2 with I7.7/I11.3. Finally, Fig. [H] shows that 
the only effect which is consistent with the observations is the variation of the fraction of 
ionized PAHs. 



1.50 



"2. 1 



0.50 



Band Ratio Calibration 
-i — 1 — 1 — I— I— I 1 1 1— 

Orion Bar 



NGC 2023 




NGC 7027 



_i I I i_ 



500 1000 3000 
(Go/nJx(T /lO^K)"-^ [cm«] 



10' 



o 

0) 

as 

1000 

(= 

o 

o 



Empirical Diagnostic 



o Dwarfs 
© Spiral 
*SB/AGN 
□ HII/PDR 
• M82 




—I I I i_ 



10 



7/1 11.; 



Fig. 9. — Empirical relation between the aromatic feature intensity ratio and the ionization-to- 
recombination ratio. Left panel: calibration of the band ratio, using well calibrated sources; the values 
of Go, and Tg^s come from detailed PDR modeling. Right panel: inversion of the right panel relation; 
the sources shown are the ones from Fig. O The grey stripe is the line fit ±1 x a. 
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We emphasize that our sample contains mainly moderately obscured star forming galax- 
ies and Galactic regions, where AGNs are under-represented. The effects t hat have bee n 



ruled out here may play a role f or peculiar s ources like elliptica l galaxies (IKaneda et al. 



2007 


), LINERs ( 


Smith et al. 


2007) 



4. THE AROMATIC BANDS AS A DIAGNOSTIC TOOL OF THE 

PHYSICAL CONDITIONS 

4.1. Empirical Calibration of the Mid-Infrared Band Ratios 

The variations of aromatic feature intensity ratios are mainly controled by the frac- 
tion of charged PAHs (Sect. 13.21) . This fraction of charged PAHs is directly related to the 
ionization-to-recombination ratio, Gg/'^e^/^gas; where Gq is the UV field density, Ue is the 
electron density, and is the gas temperature. Therefore the observed aromatic feature 
intensities can provide useful constraints on the physical conditions in the region where they 
are emitting. 
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Fig. 10. — PAH and PDR self-consistent modeling, taking into account the radiative transfer, PAH stochas- 



tic heating and charge exchange between gas and dust I^Galliano et al\\2008(\) . Left panel: SEDs at differ- 
ent locations (corresponding to different attenuations Ay) within a given slab PDR having a density of 
ngas = lO'^ cni~'^ and incident UV field of Go = 10* in Habing units. Right panel: theoretical aromatic 
feature ratio diagram; each (Gq, n^) coordinate corresponds to a single face-on integrated PDR. 
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We used the Gn 



Ue and Tgas values from the hterature of a few well-studied so urces, 



in order to provide an empirical calibration of the band ratio flGalliano et al.ll2008bl ). The 
left panel of Fig. [H] shows this calibration. It demonstrates that the I6.2/I11.3 ratio (or the 
relative fraction of PAH"*") is high for either high Gq (proportional to the ionization rate) or 
(proportional to the recombination rate). This calibration can be applied to 



low ne/i/Tj 
the sources in Fig. [3] (right panel of Fig. [9]). 



4.2. Detailed PDR Modeling of the Polycylic Aromatic Hydrocarbon Emission 



The empirical relation of Fig. [9] can be generalised theoretically, using detailed PDR 
modeling. The left panel of Fig. [10] shows the dust spectral energy distributions (SEDs) 
at different points within an homogeneous slab PDR. The radiative transfer of the incident 
ISRF on the surface of the cloud is solved, in combination with the photoioniz ation, photodis 



so elat ion and chemical reactions within the cloud, using the PDR model of iKaufman et al. 



(120061 ). In particular, this PDR model takes into account the charge transfer between gas 
and dust, and thus det e rmine s the fraction of charged PAHs at each point within the cloud. 
Using the iDraine fc Lil (120071 ) absorption efficiencies, we s olved the stochastic h eating of the 
PAHs and small grains for each Ay (left panel of Fig. [lOl iGalliano et al.l l2008d ) . 



The right panel of Fig. [TO] shows the value of the PAH band ratio integrated over each 
PDR of UV field Gq and density n^. This type of diagram is usual for PDR gas lines. 
The originality of this work is to extend it to dust features. This figure is in qualitative 
agreement with Fig. [9] It shows that the I6.2/I11.3 ratio is constant over constant Go/^e 
lines. The I6.2/I11.3 is high at the top-left corner of the panel, for high Go and low rig, where 
the PAHs are mainly positively charged. The I6.2/I11.3 is also high at the botton-right corner 
of the panel, for low Gq and high Ue, since the PAHs are mainly negatively charged. 



The present study has developed the observed mid-IR features as a quantitative tool 
to probe the physical conditions (e.g. Go/ue x ^/T^) in the emitting regions. We expect 
that this study will be of fundamental value for the interpretation of Spitzer data as well as 
future Sofia and JWST observations of galaxies in the nearby and early universe. 
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